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1  Introduction 


During  the  reporting  period  (September  1,  2010  through  January  31,  2012),  we  reviewed,  refined 
and  generalized  the  theoretical  foundations  to  analyze  and  interpret  the  optical  and  in-situ  data  that 
we  had  collected  during  a  field  experiment  conducted  in  June  2010  at  the  Boulder  Atmospheric  Ob¬ 
servatory  (BAO)  site  near  Erie,  CO.  The  optically  retrieved  time  series  of  (1)  temporal  fluctuations 
of  the  path-averaged,  vertical  temperature  gradient,  of  (2)  the  path-averaged,  beam-transverse  hor¬ 
izontal  velocity  component,  and  of  (3)  the  path-averaged  temperature  structure  parameter  turned 
out  to  agree  very  well  with  in-situ  reference  measurements.  Preliminary  results  have  been  presented 
at  various  conference  presentations  and  seminars,  and  a  masters  thesis  (Tichkule,  2011)  has  been 
prepared  on  the  basis  of  this  dataset.  Various  manuscripts  are  being  prepared  for  submission  to 
scientific  journals. 

In  May/June  2011,  a  more  sophisticated  experiment,  involving  four  sonic  towers  along  a  150  m  long 
propagation  path  was  conducted,  again  at  the  BAO  site. 

The  data-acquisition  system  and  the  spectral  analysis  routines  that  we  have  developed  over  the 
years  (supported  in  part  through  the  current  grant  and  the  previous  ARO  grant)  formed  the  basis 
of  a  pioneering  study  on  ocean-generated  infrasound  observed  at  a  location  100  miles  inland  (Sub- 
ramanian,  2009;  Subramanian  and  Muschinski,  2011).  That  study  is  only  tangentially  related  to 
the  current  research  project  and  is  not  described  here  in  any  further  detail. 

In  Section  2,  we  summarize  the  general  theoretical  approach  to  interpret  the  various  statistics 
computed  from  the  time  series  of  the  centroid  coordinates  of  the  four  test  lights  simultaneously 
observed  with  one  or  two  telescopes.  (So  far,  we  have  used  only  a  single  telescope.)  In  Section  3, 
we  present  and  discuss  two  measurement  examples:  First,  time  series  of  optically  retrieved  path 
averages  of  the  temperature  structure  parameter;  second,  time  series  of  optically  retrieved  path 
averages  of  the  beam-transverse  horizontal  velocities.  Both  time  series  are  compared  with  the 
respective  observables  extracted  from  in-situ  measurements  collected  with  sonics. 

Work  on  this  project  stopped  in  August  2012  when  the  P.I.  (Andreas  Muschinski)  left  the  University 
of  Massachusetts  and  joined  the  Boulder  office  of  Northwest  Research  Associates  (NWRA).  M.S. 
graduate  Shiril  Tichkule  followed  the  P.I.  to  Boulder  and  joined  the  Ph.D.  program  of  the  University 
of  Colorado’s  Dept,  of  Electrical,  Computer  and  Energy  Engineering  in  September  2011.  The  project 
was  re-awarded  to  the  P.I.  through  his  new  employer,  NWRA,  on  Alarch  1,  2012. 


2  Theory 


Within  the  framework  of  weak-scattering  theory  (Rytov  approximation,  Born  approximation,  or 
geometrical-optics  approximation),  the  aperture-averaged  AO  A  fluctuation  can  be  written  as  a 
volume  integral  over  a  component  of  the  refractive-index  gradient  fluctuation.  That  is,  the  vertical 
AOA  fluctuation  ocij  can  be  written  as 


aij  ~ 


dn(r) 

dz 


d3r, 


(1) 
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and  the  horizontal  AOA  fluctuation  can  be  written  as 
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(2) 


where  the  indices  i  and  j  denote  the  source  and  the  receiver,  respectively,  such  that  the  spatial 
weighting  function,  or  sampling  function,  y.y( r)  is  an  instrument  function  that  depends  on  the 
absolute  and  relative  locations  of  the  source  and  the  receiver  and  in  general  depends  also  on  other 
characteristics  such  as  the  diameter  and  shape  of  the  source  and  the  diameter  and  and  shape  of  the 
receiving  aperture.1  Here, 


r  = 


(3) 


where  the  x-axis  points  from  the  source  to  the  receiver,  the  2-axis  points  up,  and  the  y- axis  points 
to  the  left  (seen  from  the  source),  such  that  the  coordinate  system  is  right-handed. 

Therefore,  the  covariance  of,  say,  the  horizontal  AOA  fluctuation  observed  with  the  source  i  and 
the  receiver  j  and  the  mth  local  time  derivative  of  the  horizontal  AOA  fluctuation  observed  with 
the  source  k  and  the  receiver  /  is  given  by 
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or,  in  a  more  compact  notation, 
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is  the  product  of  the  two  instrument  functions  characterizing  the  source-receiver  combination  i,  j 
and  the  source-receiver  combination  k,  l ,  respectively,  and 

)(  /  //X  =  /  dn(r')  dm  dn(r")  \ 

\  dy '  dtm  dy "  j 

is  the  two-point  cross-covariance  function  of  the  y-component  of  the  gradient  of  the  refractive-index 
fluctuation  at  the  location  r'  and  the  y-component  of  the  mth  time  derivative  of  the  y-component 
of  the  gradient  of  the  refractive-index  fluctuation  at  the  location  r". 

Two-point  cross-covariance  functions  very  similar2  to  Ryy\ R^z  Ry7\ r',r"),  and 

Rzy  \ r',r")  were  introduced,  evaluated  and  discussed  by  Muschinski  (2004)  and  Muschinski  et  al. 
(2005)  in  the  context  of  a  generalized  measurement  theory  of  clear-air  Doppler  radar  windprofilers. 

1One  of  our  working  assumptions  lias  been  that  we  can  idealize  the  receiving  aperture  as  a  circular,  uniformly  filled 
aperture,  and  the  source  aperture  as  a  point  aperture.  Of  course,  these  assumptions  simplify  the  theory  considerably. 

2Muschinski  (2004)  and  Muschinski  et  al.  (2005)  considered  two-point  cross-covariance  functions  of  refractive- 
index  fluctuations  and  time  derivatives  of  refractive-index  fluctuations,  two-point  cross-covariance  functions  of  spatial 
derivatives  of  refractive-index  fluctuations  and  time  derivatives  of  spatial  derivatives  of  refractive-index  fluctuations. 
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We  have  set  up  our  experiments  such  that  the  spacing  between  the  sources  are  small  compared  to 
the  propagation  path  length  L.  Then,  using  the  Markov  approximation  along  the  propagation  paths 
and  assuming  that  the  turbulence  transverse  to  the  propagation  direction  is  locally  homogeneous 
and  isotropic  within  the  propagation  channel,  the  six-dimensional  integral  in  (4)  and  (5)  can  be 
reduced  to  a  one-dimensional  integral  along  the  propagation  direction: 


Here,  Pijki(x)  is  a  path- weighting  function  and  My™\x)  is  a  “meteorological  function”,  and  both 
functions  can  be  derived  from  Gijki{ r',r")  and  Ryy\ r',r"). 

In  the  simplest  case,  we  have  m  =  0,  and  because  in  the  inertial  subrange  Myy\x )  scales  with 
C%(  x),  the  horizontal- AO  A  covariance  tensor  can  be  written  as  a  path  integral  of  C%(x): 

( PikPji )  =  [  Pijki{x)^Cl(x)dx ,  (9) 

Jo 

where  7  is  a  constant  coefficient  that  can  be  theoretically  determined.3 

Once  Pijki(x)  and  7  have  been  derived  from  Gijki( r',  r")  and  Ryy  (V,  r"),  the  relationship  (9)  can  be 
used  to  retrieve  along  the  propagation  path  from  measurements  of  the  covariance  tensor  (fdik/Jji)- 
How  well-posed  this  inversion  problem  is  will  depend  primarily  on  the  choice  of  the  relative  locations 
of  the  sources  and  receivers.  If  one  assumes  that  C%=C%  is  constant  in  the  range  of  x  values  at 
which  Pijki(x)  has  significant  weight,  one  obtains 

{PikPjl)  =  AijklCn,  (10) 


where 

Aijki  =  /  7  Pijki(x)dx. 

Jo 


(11) 


Obviously,  this  result  can  be  used  for  an  initial  and  simple  check,  based  on  the  measured  covariance 
tensor  (/ JikPjl ),  whether  C*  is  approximately  homogeneous  along  the  propagation  path. 

The  cases  m  =  1  and  m  =  2  will  be  examined  along  the  lines  presented  in  Muschinski  (2004)  and 
Muschinski  et  al.  (2005).  We  expect  that  Myy  (x)  can  be  expressed  in  terms  of  C^(x')  and  the  lateral 
wind  velocity  vy(x),  and  that  Myy  (x)  depends  also  on  the  variance  of  the  lateral  wind  velocity. 
Correspondingly,  we  expect  Myz  (x)  to  carry  information  about  the  covariance  of  the  lateral  and 
vertical  wind  velocity  fluctuations,  that  is,  about  the  vertical  flux  of  lateral  momentum. 

The  analysis  of  the  path- weighting  functions  Pijki(x)  and  the  “meteorological  functions”  Myy\x), 
My™\x),  Mzyl\x)  and  Mz™\x)  is  to  a  large  extent  “virgin  territory”,  and  a  significant  part  of  the 
theoretical  work  will  be  devoted  to  this  analysis  as  well  as  to  the  question  of  how  stable  the  solutions 
of  the  various  inversion  problems  are. 

3Here  we  assume  that  the  diameter  of  the  telescope  is  small  compared  to  the  outer  scale  of  turbulence  and  large 
compared  to  the  inner  scale,  which  is  not  necessarily  the  case  in  the  nocturnal  ASL. 


3 


3  Observations 


Since  2005,  we  have  conducted  a  series  of  focused  field  experiments  involving  a  large-aperture 
telescope  equipped  with  a  CCD  camera  and  up  to  twelve  sonics  mounted  on  up  to  four  portable 
towers  along  the  propagation  path  (Behn,  2006;  Cheon  et  ah,  2007;  Cheon,  2008;  Behn  et  ah,  2008; 
Whiteman  et  al.,  2008;  Tic.hkule,  2011),  and  up  to  three  ultrasensitive  quartz-crystal  barometers 
in  order  to  observe  the  mean  pressure  and  pressure  fluctuations  associated  with  gravity  waves  and 
infrasound  (Subramanian,  2009;  Subramanian  and  Muschinski,  2011). 

In  the  following  two  subsections,  we  present  measurement  examples  in  order  to  demonstrate  that 
temperature  structure  parameters  and  lateral  wind  velocities  can  be  effectively  retrieved  from  optical 
AOA  statistics  estimated  from  test-light  image  centroids.  We  will  show  that  the  optically  retrieved 
quantities  (based  on  the  geometrical-optics  approximation)  compare  well  with  independent  in-situ 
observations  even  if  the  data  are  processed  in  a  relatively  simple  manner.  These  measurement  exam¬ 
ples  are  applications  only  of  simple  cases  (single  telescope,  m=  0,  m=  1)  of  the  general  measurement 
theory  described  in  Section  2. 

We  collected  the  data  presented  in  the  following  in  a  field  experiment  conducted  during  the  night 
of  15/16  June,  2010  at  the  BAO  site.  The  experimental  setup  was  practically  identical  with  the 
one  described  by  Cheon  et  al.  (2007),  except  that  LEDs  were  now  used  instead  of  krypton  flashlight 
bulbs,  and  twelve  sonics  (instead  of  a  single  sonics)  were  mounted  on  two  20-ft  towers  along  the 
propagation  path. 

The  propagation  path  was  horizontal,  at  1.7  m  AGL,  and  182  m  long,  with  the  telescope  at  one  end 
and  the  test-light  array  at  the  other  end.  A  CCD  camera  captured  telescope  images  of  four  LED 
test  lights  (vertical  spacing  7.5  cm,  horizontal  spacing  10  cm)  at  a  rate  of  30  frames  per  second. 
From  each  image,  vertical  and  horizontal  centroid  coordinates  were  estimated  for  each  of  the  four 
test  lights,  such  that  each  image  provided  four  vertical  (cci, 0:4)  and  four  horizontal  (/?i, ...,  ^4) 
angles-of-arrival  (AO As). 

Six  sonics  were  mounted  on  a  20-ft  tower  located  47  m  away  from  the  telescope  and  2  m  away  from 
the  propagation  path.  The  two  lowest  sonics — the  only  ones  that  are  used  for  the  measurement 
examples  presented  here — were  at  1.45  m  and  2.13  m  AGL.  The  sampling  frequency  of  the  sonic- 
data  (temperature  and  the  three  wind-velocity  components)  was  32  Hz. 


3.1  First  measurement  example:  Temperature  structure  parameter 


Figure  1  shows  two  time  series  of  the  temperature  structure  parameter,  Cj.,  one  of  which  was 
estimated  from  sonic  data  (black  dots)  and  the  other  one  from  horizontal- AO  A  variances  (red  dots). 
Each  data  point  is  a  30-s  estimate. 

The  optically  retrieved  Cj,  values  were  obtained  by  means  of 


C\  = 


(ft  -fhypl/3 

7 sbL 


(12) 


where  the  overbar  stands  for  averaging  over  30  s,  D  =  35  cm  is  the  aperture  diameter,  js  = 
(3/8)  x  2.383  =  1.064  is  the  coefficient  for  spherical  waves  in  the  geometrical-optics  limit  (Cheon  and 
Muschinski,  2007,  eq.  37  on  p.  420),  L  =  182  m  was  the  propagation  path  length,  and  b  =  ( ap/T 2)2 
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is  the  constant  of  proportionality  that  relates  the  temperature  structure  parameter  to  the  refractive- 
index  structure  parameter,  where  a  =  7.9  X  10-7K/Pa  for  red  light,  p  =  839  hPa  was  the  surface 
air  pressure,  and  T  =  288  K  was  the  air  temperature. 

The  sonic-retrieved  Cj  values  were  obtained  by  means  of 


(ap\  2  (Ari2-ATi2)2 

VTV  A  2/3 

^12 


(13) 


where  AT12  =  T{z2)  —  T{z\)  and  Azu  =  Z2  —  z\  and  where  z\  =  1.45  m  AGL  and  Z2  =  2.13  m  AGL 
were  the  heights  of  the  lowest  two  sonics.  (The  height  of  the  propagation  path  was  1.7  m  AGL.) 


Mountain  Daylight  Saving  Time  MDT  [hh:mm] 


Figure  1:  Time  series  of  the  temperature  structure  parameter,  Cj,  1.7  m  AGL,  measured  with  a 
pair  of  vertically  spaced  sonics  (black)  and  retrieved  from  horizontal  AO  A  fluctuations  of  a  single 
testlight  observed  with  a  CCD  camera  mounted  on  a  telescope  at  a  distance  of  182  m. 

Figure  1  shows  that  the  two  Cj  estimates  track  each  other  nicely  down  to  time  scales  of  30  s, 
even  though  the  optically  retrieved  C\  values  are  path  averages  while  the  sonic-measured  Cj,  values 
are  essentially  point  measurements.  The  sonic-measured  values  are  systematically  larger  than  the 
optical  Cj,  retrievals,  by  20  or  30%.  This  bias  may  be  caused  by  various  mechanisms,  including 
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(1)  anisotropy  of  the  turbulent  temperature  field  in  the  nocturnal  ASL  at  the  sonic  separation 
length  scale  z  12  =  68  cm,  (2)  the  fact  that  close  to  the  ground  decreases  strongly  with  increasing 
height,  and  (3)  possible  systematic  inhomogeneities  of  along  the  propagation  path  due  to  possible 
roughness  inhomogeneities. 


3.2  Second  measurement  example:  Horizontal  wind  velocity 

Figure  2  shows  time  series  of  beam-transverse,  horizontal  wind  velocities  observed  with  a  sonic 
anemometer  (black  dots)  and  retrieved  from  the  delays  of  horizontal- AO  A  fluctuations  from  a  pair  of 
horizontally  spaced  test  lights  (LEDs)  observed  with  the  same  telescope  (blue  dots).  The  horizontal 
spacing  was  10  cm. 
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Figure  2:  Time  series  of  path-transverse,  horizontal  wind  velocities  observed  with  an  ultrasonic 
anemometer  (black  dots)  and  retrieved  from  optical  AOA  fluctuations  by  means  of  the  slope-at-zero- 
lag  method.  Each  data  point  represents  a  10-s  average. 

The  transverse  wind  velocities,  vt,  were  retrieved  by  means  of  the  “slope-at-zero-lag”  method.  Details 
of  the  data  processing  can  be  found  in  Shiril  Tichkule’s  MS  thesis  (Tichkule,  2011). 

Each  data  point  represents  an  average  over  10  s.  The  two  time  series  track  each  other  well  but  not 


6 


perfectly,  which  is  not  surprising  because  the  optical  data  represent  (weighted)  quasi-instantaneous 
path  averages  of  the  transverse  wind  velocity,  while  the  sonic  data  are  quasi-instantaneous  point 
measurements.  The  r.m.s.  difference  between  the  two  time  series  is  about  20  cms"1,  which  is 
probably  the  result  of  spatial  variability  in  the  horizontal  wind  field  at  length  scales  between  a 
few  tens  of  meters  and  the  path  length.  However,  there  appears  to  be  no  significant  bias  between 
the  optically  retrieved  and  the  sonic-measured  velocities.  2-h  averages  of  the  differences  between 
sonic-measured  and  optically  retrieved  wind  velocities  obtained  in  this  experiment  do  not  exceed  3 
cms-1. 

4  Summary  of  the  most  important  results 

We  demonstrated  that  the  optical  angle-of-arrival  fluctuations  observed  with  a  large-aperture  tele¬ 
scope  equipped  with  a  moderately  fast  CCD  camera  can  be  used  to  retrieve,  with  high  sensitivity 
and  accuracy,  path  averages  of 

•  the  vertical  temperature  gradient  (not  discussed  in  this  report) 

•  the  transverse  wind  velocity 

•  the  temperature  structure  parameter 

along  a  horizontal  propagation  path  (200  m  path  length,  1.7  m  above  ground  level)  in  the  atmo¬ 
spheric  surface  layer. 
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